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ABSTRACT 

In tile present work, a methodology of interac- 
tive, design and analysis of multi-speed machine tool gear 
box shafts has been developed, Ctomputer program developed 
here gives an appropriate design of gear box shafts 
selected from the different specified options. The 
program can also analyse general machine shafts, when 
geometric properties, material properties and general 
loadings are given as inputs. This analysis procedure 
provides complete information about the deflections, 
stresses, and the factor of safety with considerations 
of stress concentrations, failure due to yielding and 
fatigue. The program includes graphics display of 
results so as to help the engineer in the process of 
evaluating the design. The analysis procedure is based 
on transfer matrix method using the fundamental equations 
of motion for the bending of a beam. The proposed 
approach has been illustrated vrith the help of four case 
studies and ihe results obtained have been discussed. 



CHAPTER - I 


INTRODUCTION 

1 • 1 RQle of Gear Boxes : 

The gear boxes of general purpose mac±iine tools 
provide a wide range of cutting speeds and torques from 
a constant speed power input enabling proper cuttirjg 
speeds or torques to be obtained at the spindle as 
required in the case of spindle drives and the desired 
feed rates in the case of feed drives. Optimum cutting 
speeds and feed rates enable the operator to obtain 
optimum rate of metal ranoval and minimum operation time. 

Witii a constant speed power source there is a 
need for some method of varying the speed over a given 
range. Stepless mechanical and electrical drive can 
provide infinite speed variation. However# the torque/ 
speed characteristics of available stepless drives do 
not meet the requirements of spindle drives which demand 
an increased driving torque at lower output speeds in 
order to maintain a constant rate of metal removal. The 
stepless drives which do possess the required torque/ 
speed characteristics are limited by the speed range 
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over whicih these characteristics can be maintained. In 
order to provide a wide range of opera tiing speeds together 
with adequate toirgue at lower spindle speeds# it is 
necessary to use gear boxes which enable the required 
spindle speed range to be covered in a number of dis- 
crete steps. 

1.2 Design of Shafts in Gear Boxes : 

The gear box consists of sliding gears# non- 
sliding gears# shafts# bearirjgs and casing. Thus shafts 
are the integral part of a gear box and the design of 
these shafts is of prime importance. 

Shafts are rotating members whose function is 
to transmit power and motion. Since shafts are so 
commonly used, it is very important to employ a good 
design procedure. Shaft design consists primarily of 
determining the correct shaft diameters to provxde satis- 
factory strength and rigidity. Most shafts are gene- 
rally subjected to fluctuating loads of combined bending 
and torsion with various degrees of stress Concentra- 
tion. 

Since shafts are the most commonly designed 
machine elements, the assumptions and sequence of design 
steps are instrumental to the results of the final design. 

must tahe into account applied forces# 


This design process 
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permissible deflection limits and stress concentrations 
at critical sections and locations such as keyways, etc 
Most of the shafts can be treated as beams which enables 
usage of strength of material equations for deflection, 
slope, bending moments, and shear force along the shafts 

[1,2], 

The gears must be suitably mounted on shafts 
of such size as properly to take care of maximum torque 
loads imposed by rated hors^ower and speeds. Theore- 
tically these shafts should be capable of withstanding 
as much load as can be transmitted by the gear teeth and 
must also be rigid enough to resist any bending loads 
that might possibly tend to affect proper engagement of 
the gear teeth. The shafts must also be provided with 
the necessary keyway s, splines- . etc. by means of which 
the motor and driven units can be connected. 

Since horsepower is a function of the piroduct 
of speed and torque, the high speed shaft is small and 
the low speed shaft is larger so as to give an increase 
in strength which is in proportion to the decrease in 
speed. 

The shafts must be supported on bearings, usu- 
ally of the ball or roller type, of such size, type and 
capacity as properly to care for the loads, both radial 
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and thrust, which would be setup under maximum operatinor 
conditions. 

^ ^ Qb-lective of the Present Work ; 

With computers the design process of shafts 
or any ma<±iine element can be significantly improved. 

To be able to effectively utilize the computer, a correct 
and simple algorithm for shaft design must be developed 
in detaij. By employing the computer, the results of 
minor adjustments or major changes in a design can be 
quickly assessed. This helps the engineer to test scy 
different designs without doing large amount of repeti - 
tive calculations manually. Also using interactive 
graphics, the engineer can examine the possible models 
on a computer terminal before executing the solution 
procedure. 

Gear boxes with a large nxamber of steps within 
a given range would be bulky and expensive. Hence thes^ 
should be so designed that while fulfilling the func- 
tional requirements, the assembly should be economical 
to manufacture. The cost of the gear box is related to 
the number of shafts and bearings required and the total 
number and size of the gears. From different possible 
arrangements of the gears, the layout which gives com- 
pact size and lower cost while still fulfilling the 
technical requirement of the system should be chosen 

[3,4,5]. 



5 


To reduce the cost of the gear box the optimsT 
sizes of the shafts should be founa out from different 
possible arrangements of the gears. 

Once the number of gears, shafts and the rambe' 
of teeth on each gear and also the input horsepower are 
■known, one can searc±i for the possible arrangements of 
the gears. In order to get minimum length and diameter 
of a shaft it should be kept in mind that the distance 
between two gears of the fixed block must be equal to 
two gear widths. Thus by giving possible different go^ 
arrangements the loads on the intermediate nonsliding 
shafts can be calculated. By knowing the loadings, 
geometric properties such as diameters and bearing j^c 
tions and material paxjperties, the factor of safety' 
deflection along the shaf’'- can be deterrr ned. This is 
essentially an analysis problem which determines how Su-x. 
a shaft is for the proposed loading. Afuer analysis c 
each possible gear arrangement, the best possible design 
of a shaft is determined. 

The aim of this ■'jork is to develop the solution 
procedures and algorithms for the suggested problems. 
Since stress concentrations cannot be avoided, the algo- 
rithms incorporate stress concentration to provide an 
accurate solution technique. The major aim of this work 
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is to supply a reliable solution procedure. This per- 
mits an engineer to devote more time to the engineering 
aspect of shaft design and less time on repetitive cal- 
culations. 

1,4 Scope and Limitations : 

Most shafts are statically indeterminate systems 
which require both the equations of statics and deforma- 
tion to be solved. These systems can be analysed by 
using finite difference/ finite element/ or transfer 
matrix methods. Using any of these methods enables 
modelii^ of shafts with complex loading and arbitrary 
cross-sections. But the present work deals with only 
circular cross-sections. 

The transfer matrix method is chosen over the 
other two methods for 'several reasons [f]* The finite 
element method gives an accurate solution/ but its solu- 
pj-Qcgdure is technically more complex. Concerning 
programming capabilities/ the finite element method 
requires greater amount of computer storage and is itiore 
difficult to program. When programming the finite diffe- 
rence method/ more storage and more computer time is 
needed than for transfer matrix method. The transfer 
matrix method provides an efficient procedure for find- 
ing the deflection, slope, bending moment and shear 
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force at any location along the shaft [6,7], A major 
advantage of this method is to model all possible types 
of loads, supports, and geometry, with a simple solution 
procedure that can be easily programmed. For exairple, 
a stepped shaft with different material properties and 
loadings along the shaft can be easily analysed. The 
possible types of loadings considered in this work are 
concentrated loads, distributed loads, torques, concen- 
trated moments and axial forces. Supports can be located 
anywhere along the shaft with no theoretical limit on 
the number of supports. Rigid supports, translational 
and rotational springs also can be modeled. This peinmits 
rigid, pinned, flexible, and free boundary conditions to 
be considered. 

The major types of stresses to which the shaft 
is subjected are shear stresses caused by torsional load- 
ing and shear forces, and normal stresses caused by axial 
loading and bending [l,S,lo]. The shear forces and 
bending moments are computed by transfer matrix analysis, 
where as the torques and axial loads are known input 
data. The state of stress at each point along the shaft 
must be combined to form an equivalent stress by employ- 
ing an appropriate failure theory. The types of failure 
theories commonly considered in design texts are the 
maximum normal stress theory, maximum shear stress theory 
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and -the distorsion energy theory [s,9,10]. The maximum 
normal stress theory is generally recommended for brittle 
material and is not usually considered for shafting* So 
only the last two theories are considered here. 

Since some shafts are subjected to fluctuating 
loads with nonzero means, the resulting stresses are 
also fluctuating with nonzero mean stress. However, most 
of the fatigue data available is from tests with zero 
mean stress i.e, completely reversed load. Therefore, 
it is important for the design equations to include a 
stress which has an equivalent zero mean cyclic stress. 
Several methods have been devised for estimating the 
equivalent zero mean stress when given general stress 
condition. These methods include the modified Goodman, 
Soderberg and Gerber relationships [9]. The Soderberg 
relationship is a conservative approximation and is one 
of the most commonly applied methods. Thus it has been 
incorporated in this work along with the two failure 
theories to calculate an equivalent or resultant uni- 
axial stress. Once the equivalent stress is found, the 
failure stress can be divided by the equivalent stress 
to produce the factor of safety at each desired loca- 
tion. Prom the factor of safety at the different stress 
concentrations and at other critical points along the 
shaft, the smallest factor of safety for a given shaft 
can be determined. 
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Shaft design requires determining the shaft 
diameters that ensure a safe design. Thus knowing all 
the loading, support conditions and the material pro- 
perties the diameter of the shaft can be calculated. 

Thus in case of gear box the diameters and the length? 
of the shaft for each possible different cases are found 
out and the smallest and best of all is chosen as idle 
optimal shaft size. 



CHAPTER «• 2- 


TRANSPER MATRIX METHOD 

2.1 I ntro ducblo n : 

For the design of rotating shafts# the shaft' s 
strength and rigidity are important. Strength evalua- 
tion requires the calculation o£ stresses and stiffness 
considerations require the calculation of slope and 
deflection. 

To calculate stress along the shaft the bending 
moment and shear force must be known. Stiffness of shaft 
is of interest because the slope of shaft can be a limit- 
ing case near the bearing support and excessive deflec- 
tion might cause improper engagement of gears and inter- 
ference with another machine parts. 

By transfer matrix we can calculate simultaneously 
the deflection# slope, bending moment and shear force at 
any point along the shaft. In this work the differential 
equations developed assume small deflections and linear 
elastic material properties. 

The transfer matrix approach for modelling of 
rotating shafts and beams as given by Pilkey and ChangE?] 
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is presented in this c±iapter. Thus the formulation deals 
with static stability analysis of beams by considering 
arbitrary loadings and several types of supports. 

2.2 Governing Equations for Shaft Analysis : 

Since a shaft can be assumed to be a rotating 
beam, the equations of motion for a beam as presented 
by Timoshenko [l] can be used to develop the transfer 
matrices. These general equations for the bending of 
beams as given by Pilkey and Chang [?] are. 
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In this work presented here, these equations can be 
simplified by number of assumptions such as, 

i) The shaft rests on a rigid foundation. Thus 
K* = = O 

ii) Since shafts are rotating, a radial thermal 
moment is unlikely to develop. Thus « O 
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iii) Since a moment, intensity or distributed moment 
is not a common loading and which is not con- 
sidered in this work. Thus C (x,t) = 0 

iv) The axial loading effect on the motion of a 

beam is not significant for small deflections 
of shafts p Thus P = O 

v) When considering a Euler Bernoulli beam i.e* 

uniform ordinary beam the rotary inertia 

2 

term is neglected. Thus r -—r- = O 

vi) Since a static analysis is being Gonsidelred, 
time dependence is also eleminated. 

Considering these assumptions ,the above ogtia- 
tions 2,1 to 2.4 become 
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where, A 
K 


Cross sectional area 
Shear form factor (SFP) 


This shear form factor K for some standard shapes can 

s 

be found in reference [ 11,12]. For circular cross 
section it is 

£1 / -f » ^ 

( 2 . 10 ) 
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(i + y) 

7“+ 6 r 


where, y 


Poisson's ratio 


After rearranging and integrating. Equations 2.5 to 2.8 
become 
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The loading functions, , P^ are given for 

general types of loading in Tabl e 2.1 . These loading 
functions are for the general solution of a shaft 
that is supported on the ends. If more than one loading 
occurs on the shaft, the loading function for each type 
of loading is added using superposition. The table uses 
singularity functions whic±t state that for n > 0 

0 if X < a 
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Equations 2,11 to 2.14 can be written in matrix notation 
as. 
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or in extended form as 
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Table 2.1 for loading functions , Pj^ and F^ for 

general uniform shafts a^s shown in table 2,1 . The column 

matrices are known as state vectors, since they contain 

the state variables y, 9 , M and V, The state vector contain- 

ino V / 9 1 M and V is located at the left end of 

^-^o o o o 

the shaft. The square matrix is called the transfer 
matrix, since it 'transfers' the state vector at the 
left along the shaft. This transfer matrix can be re- 
presented in general form as. 
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where each term in [Uj_] gives the influence of the indi- 
cated left end state variable on the indicated right end 

state variable. For example, U gives the influence 

yw 

of 9 ^ on y and is equal to - x in Equation (2.17). 

When the shaft has a complicated loading or 
in-span supports, the simultaneous equations become 
complex and difficult to manipulate into a specific 
solution. However, the solution process for complicated 
shafts can be easily solved with the transfer matrix 
method . 

2.3 Sign Convention ; 

Special attention must be given to the sign 
convention used in the development of the different types 
of transfer matrices, to give the solution. Using a 
right handed coordinate system the x axis is positive 
to the right, z axis is positive out of the paper and 
they axis is positive downward as in Figure 2.1. All 
external forces are positive if applied in the direction 
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of a positive axis* For example, a concentrated force 
applied in the direction is positive downward. Like- 
wise, internal forces and moments acting on a positive 
face (face on the right end of a free body diagram, see 
Figure (2.2) )are positive when their vectors are in the 
direction of a positive coordinate* Internal forces and 
moments acting on a negative face are positive in the 
direction opposite the positive coordinate, as shown in 
Figure (2.2)* Following this right hand rule, the def- 
lection and slope are also positive if their vectors are 
in the positive coordinate direction as in Figure (2*1). 
2.4 Types of Transfer Matrices : 

There are two types of transfer matrices, field 
matrices and point matrices. 

i) Field Matrices: These matrices are generated for 

finite lengths along the shaft wiidi constant material 
properties and constant geometrical properties. The field 
matrix^ for section of length 1 is given by as, 

1 -1 

0 1 

0 0 

^ i 0 0 

0 0 

i Derivation of fleld ' "m atrix is given in Appendix I. 
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Fig, 2.2 Sign convention for positive bending. 
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The general form of a distributed load is given in 
Figure(2,3). The loading equations can be obtained 
from table 2,1 by letting x t= a 2 and 1 = a 2 - and 
superposing the cases for a constant distributed load 
and a triangular distributed load. These equations for 
the loading functions reduce to 
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( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


where w 

AW 

b-1 


r > 

s= constant distributed load 
= gradient of distributed load (/uw is the 
change in the loading and ts 1 is the change 
in length) 

If section is not loaded, the loading functions are set 
equal to zero. 

Field Matrices for Axial Load ; These matrices are also 
generated for finite lengtli along the shaft with constant 
material properties and constant geometric properties. 
Field matrix for any section of length 1: for axial load- 
ing in given by [u^l as. 




Fig 2.3 Simply supported shaft with linearly 
varying distributed load. 
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P axicl force is taken as the actual axial force 
in the segment as found fix)m axial equilibrivim 
requirements . 

P-P^tt Matrices ; These matrices take into account 

the concentrates occurrences. Here a point occurrence is 

defined as a concentrated load or moment. The cases of 

concentrated occurrences considered in this work are 

point loads, concentrated moments, in-span supports, 

translational springs and rotational springs (and lumped 

masses). The point matrix[,U^], is presented for all 

of these cases except for in-span supports in Table 2.2. 

Noticing that for each type of point occurrence, there 

is only one entry to the identity matrix. For example 

the point matrix for a translational spring is formed 

by substituting the value of the spring constant for 

the (4,1) element in the (5x5) identity matrix given 

in Table 2 *2 and where , c and W are zero. Since the 
2 

term Ip (v-f in Table 2.2 is small compared to oihers, 
thus in the program developed ^this term is not included. 
Usually the shafts are supported on bearings and the 
boundary condition is given as simply supported one. 

But if we know the mass of the bearing pedestal, stiff- 
ness coefficients and damping coefficients and whirl 
frequency then the boundasry condition at the bearing can 



POINT MATRICES 
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be given as free one but there will be a point matrix 


[ 7 , 14 ] which is as shown. 



where 

(K. + sc. ) (K. + SC^ + M. ) 

r, 11,2 2 1 

~ ^ 

Kj + SC^ + ^2 + SC 2 + 


= Stiffness coefficient for bearing (forc^lengtl 
K 2 =s Stiffness coefficient for bearing Cforce/length ) 

= Damping coefficient for bearing (force time/lengi' 

C 2 SB Damping coefficient for bearing (force time/1 eng tli'5 
= Mass of the bearing pedestal 
S = i Q 

Q = VThirl frequency (radians/time) 


2 .5 Global Matrix; 

After the transfer matrices are formed, the global 
matrix is formed by multiplying together the individual 
transfer matrices. The global matrix can be defined as, 

which after matrix multiplication becomes 


(2.27) 
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where [u] is a global matrix. Here/ is the 

state vector which contains V , © / M and V . For 

^ o o o o 

the case in Equation (2,26)/ [u^] is the first field 
matrix at the left end of the shaft and is a point 

matrix that occurs at the right end of field matrix [u^]. 
The transfer matrices in Equation (2.26) are so arranged 
that first matrix is the right most matrix . Then 

matrix which is just to its left and so on till the 
left most Aatrix is reached. The matrices must be multi- 
plied together from left to right, since matrix multipli- 
cation is not commutative. For example, when multiplying 
[u^] the point matrix must be on the left side to 

obtain the correct answer. Once the global matrix has 
been found, the boundary conditions can be applied. 

2,6 Boundary Conditions ; 

In this worb three types of boundary conditions 
considered are fixed ends. Pinned ends and free ends. 
Pinned and fixed ends must have a rigid supports located 
at the specified end. A fixed end is a condition for 
whidh the deflection and slope are equal to zero. In 
case of pinned end the deflection equal to zero, but 
since a pinned end is free to rotate, the bending moment 
is equal to zero at that location. In case of free 
end, both the bending moment and shear force are zero. 

If a shaft has a spring support or point load at the 
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end of an overhung section, this occurrence is modelled 
as a point matrix on a free end* 

The purpose of applying the boundary conditions 
is to solve the governing equations and to obtain the 
state Vector at the left end* Two state variables are 
always known at the left end from the boundary conditions, 
and the other two state variables are calculated from the 
boundary conditions at the right end. Two equations with 
the unknown state variables are formed from the global 
matrix and the four boundary conditions. A simple pro- 
cedTire to calculate these equations has been developed 
by Pilkey [v], and the same is used in this work. In 
this procedure, the columns of the global matrix are 
cancelled for the state variables that are zero at the 
left end. For example, if the deflection is zero at the 
left end, the first column of the global matrix is can- 
celled. Rows are cancelled where the state variables 
at the right end are unknown. After this two equations 
remain which contains two unknown state variables at the 
left end can be solved and values of the unknowns can 
be found. * 

To illustrate this procedure, a shaft is con- 
sidered with a pinned end at the right end and with a 
fixed end at the left end, as shown in Figure (2.4). 

X A problan is" solved in Appendix II, to illustrate 
the procedure- 





Fig 2.4 Loaded shaft used for application 
of boundary conditions. 
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The 'transfer matrix notations for this shaft becomes^ 

{®|x=a^ = [^ 3 ] tu,] [Uj] [U^] [Uj] (s|^^ (2.28) 

The general global matrix obtained after matrices multi- 
plied together in Equation (2.23) is as. 


Y = 0 
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M = 0 


V 


1 ■ 

X=ij 


^YY .He ^vM ^yV 

! ■ * 

^MO 

■’^y" "^ve ■ ^VM ■'^W "“^y 


= O 

o 


4 q 


1 x=o 

(2.29) 


Since ■the deflection and slope are zero at the left endj 
columns 1 and 2 of -the global matrix are eliminated. 
Similarly# wi-th and unknowns at the right 

end# where L is -the length of the shaft, rows 2 and 4 
of the global matrix are eliminated. This leaves the 
following two equations 

° “ V“° V ''o 


° = + “mV '^o + 


(2.31) 


Solving these two equations for and V^# yields the 
unknowns in the state vector at the left end. The 
initial parameters i.e. the unknowns state vector para- 
meters can be calculated directly from Table 2,3 with 
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the help of boundary conditions for shafts without in- 
span supports. 

2 *7 Calculation of Deflection, Slope/ Bending 
Moment and Shear Force Along the ^aft s 

Once the state vector, jf has been found, 

the deflection, slope, bending moment and shear force 
can be calculated anywhere along the shaft. Between 
sections, results are computed by appropriately adjust- 
ing the lengti variable in the field matrix for the 
point of interest. To calculate the deflection, slope, 
bending moment and shear force at a point of interest, 
the transfer matrices from the left end Pf the shaft 
to the point of interest are multiplied together with 
the state vector at the left end. For example in 
Figure 2,4./ the state vector at x = a 2 is equal to 

{=] X=a3 = '2.32) 

While programming this procedure, it is not 
necessary to store all of the matrices. When initially 
calculating the transfer matrices, these matrices are 
multiplied together as they are developed. Significant 
parameters, such as moment of inertia, and modulus of 
elasticity, are stored for the later regeneration of the 
field matrices. No extra parameters are stored for -the 
point matrices, since all of -the required information is 
input da-ta. 
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Once the global matrix is developed the state 
vector at left end of the shaft is calculated according 
to procedure developed using boundary conditions# I'ilhen 
making a second pass to compute the forces and deflec- 
tions along the shaft/ the field matrices are regenerated 
to the appropriate length to determine the deflection, 
slope, bending moment, and shear force at the desired 
locations. However, when encountering a new field matrix, 
the state vector at the beginning of the new field matrix 
is calculated and is stored. This new state vector is 
used to calculate the deflection, slope, bending moment 
and shear force from that point on the shaft to the next 
new field matrix [l43* For example in Figure 2.4# the 
state vector at x = a 2 </and x = a^. m equation form 
this can be written as 

where [Upl is the third field matrix having an appro- 
priate adjusted length. This process saves remulti- 
plying the matrices for the second sweep across tne 

shaft. 

2.8 In-soan Supports ; 

A rigid support which does not occur at the 
end of the shaft is called an'ln~spanin3eterminate con- 
dition, The major differences in solving a shaft wrth 

/ is used to calculate ihe deflection, slope, bending 
^ moment and fdrce between x = aj 
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an in— span support are the steps in the development of 
the global matrix and the application of the boundary 
conditions. 

The solution procedure for in-span rigid supports 
is somewhat parallel to the case without in— span rigid 
supports. First, all the transfer matrices are calcu- 
lated except for the in-span point matrices. Then, the 
subglobal matrices are determined from the calculated 
transfer matrices. Each subglobal matrix contains all 
of the transfer matrices between the left end and the 
first in-span support, between each in-span support, 
and between the last in-span support and the right end. 
Thus the number of subglobal matrices is equal to one 
more than the number of in-span supports. 

For example, if two in— span supports are present 
along the shaft, then there are three subglobal matrices. 
The first subglobal matrix contains the transfer matrices 
from the left end to the first in-span support at x = L^, 
as shown in Figure (2.5). The transfer matrices are 
multiplied together using the same method as for the 
global matrix. The second subglobal matrix contains the 
matrices from just to left of the first rigid xn-span 
support bo the second in-span support at x « L^. The 
final subglobal matrix is determined in the same manner 






Fig 2.5 Shaft used for determination of 
subglobal matrices. 
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from X = L 2 to the right end of the shaft. After the 
subglobal matrices have been calculated, the point 
matrices for the in-span supports can be determined. 
Solving for the point matrix at the in-span support is 
essentially equivalent to evaluating the reaction force 
at that location. The general form of the point matrix 
for the in-span supports is shown in table 2.4, To 
solve for this point matrix, the procedure uses a boun- 
dary condition to the right of the in— span support. If 
another in-span support occurs to the right of a given 
in— span support, the bnown condition is the deflection, 
and the subscript in Table 2.4 is equal to y. However, 
if no in- span support occurs to the right of an in- span 
support, then subscript in Table 2.4 is equal to one 
of the boundary conditions at the right end. The terms 
in the point matrix for in-span support are taken from 
the subglobal matrix to the right of the in-span support. 
The global matrix is equal to 

where [ and [u^] are point matrices for in-span 

supports. The subglobal matrices are 


Subglobal matrix ■ 

1 

; Transfer matrices 

[Uj] 

2 

[u,] [ a,] [U3] [U3] 

3 

t 05] K 3 
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TABLE 2.4 

Point raatrix for in-span rigid support 

Segment i Segment i + 1 

~ 2 . 

J 



Definition for rigid support: 


s = y if there is a rigid in-span support to the right 

of X = a£. if there no in-span support to the right 

of a. . then e is one of the state variables (y^, © , M , 
i ^ ^ 

or V) that is zero at the right end of the beam i.e. to 
say that Sja is one of the boundary conditions at right 
end of the beam. are the sutglobal matrix elements 

for the segment between x » and the next in-span 
support or the right end. 
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Once the in-span point matrices are found, the 
global matrix is calculated by matrix multiplication 
in the order as discussed in Section 2,5. Since one of 
the boundary conditions at the right end is used for 
calculating the point matrix at the last in-span indeter- 
minate condition, the method of applying boundary condi- 
tion in this case is somewhat different from without 
in-span support. The four boundary conditions needed 
to solve the system eguations are two at the left end, 
one at right end unused boundary condition and one at 
the first in-span support. The two simultaneous egua- 
tions for calculating the two unknown parameters at the 
left end are formed in similar way that of in Section 
2.6. One equation is developed from the global matrix 
and the unused boundary condition at the right end. The 
other equation can be formed from the first subglobal 
ni 3 ^-trix and the known deflection at the first in— span 
support. After the state vector at left end is found, 
determining°the deflection, slope, bending moment and 
shear force is done usirg same procedure discussed in 
Section 2*7. 

2,9 Tran sf er Matrices for Torsio n 

Torsion is the result of one or more forces 
twisting a bar about a longitudinal axis. The net 
result of twisting loads, acting on a metiber is a couple 
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which has the dimensions of force and distance# The 
force unit is usually stated first to distinguish a torque 
from a bendirg moment which has the same dimensions. 


If a transverse force is applied at such a 
point in the cross-section that twistirg results^ two 
equal and opposit forces act at the center of twist, 
with the result that a bending force and a couple act 
on the section. The center of twist is located at the 
shear center when a transverse force causes twisting. 

For analysis of gear box shafts the torsion 
plays an important role. 


The general fom of transfer matrix for torsion is as 



T<^ 



U, 


TT 


0 



( 2 . 34 ) 


The field matrix for torsion is as 



1 __ 

GJ 

1 

0 



( 2 . 35 ) 


and Fm the loading functions are zero since 

^ i 

only concentrated occurrence of torque is assxamed, which 
is taken care in point matrix. 
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The point matrix for torsion is as 



1 

j-p 

0 


0 

1 

0 


0 

-T. 

X 

1 


(2.36) 



Concentrated applied torque 



Disk occurrence 



2 


The boundary conditions in this case are only two^one 
is fixed and free. Since pinned end is free to rotate 
the boundairy condition in case of pinned end is given 
as free end for torsional calculation. At fixed end 
angle of twist is zero and at free end torque is zero. 
The global matrix is found in similar way using the 
procedure discussed in Sections 2.5 to 2.7, And by 
applying boundary conditions the state vector £ S 
is found out. 


1 


x=o 



(2.37) 


Then the angle at twist and torque anywhere along the 
shaft are found out using same procedure developed in 
Section s 2.5 to 2.7. 



CHAPTER - 3 


STRESS ANALYSIS 


3.1 Introduction; 


The intention of this chapter is to introduce 
the guide lines for finding the stresses on a rotating 
shaft and thus the least factor of safety. In consider- 
ing the three dimensional loading on a shaft# the forms 
of stresses generated are shear stresses and normal 
stresses. The loadings considered on the shaft in this 
present work are^radial loading in two planes# axial 
loading and torsion. The shear stresses are developed 
by the torsional loading and radial shear forces (trans- 
verse shear). Likewise# the normal stresses are generated 
by the axial loading and the bending moments,. 


A mechanical element (say shaft) may be stressed 
in tension> compression or shear# either singly or in 
combination, by loads which are either static, fluctuat- 
ing or composed of static and variable components. Depen- 
ding upon the nature and magnitude of the loads# failure 
may occur by yielding# cracking or rupture. Statically 
loaded structures will normally perform satisfactorily 
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until the yield point of the material is reached and in 
some cases, until rupture occurs. Predictirg failure 
for a shaft or any other mechanical element with uniaxial 
loading is relatively easy, However^for the possible 
loading on shafts, failure prediction is sufficiently 
more difficult because of a multiaxial state of stress. 
Several failure theories have been advanced for predict- 
ing failure under biaxial -and triaxial stress conditions 
f 3 Com uniaxial material properties [ll]- Some assumptions 
are made while calculating stresses, bending moments and 
shear forces. 

1) The beam (shaft) is of homogeneous material 
which has the same modulus of elasticity in 
tension and compression, 

2) The shaft is straight. 

3) The cross section is uniform. 

4) All loads and reactions are perpendicular to 
the axis of the beam and lie in the longitaoinal 
plane of symmetry. 

5) The maximum stress does not exceed the prepor- 
tional limit. 

The beam is loaded only by equal and opposite 
twisting couples, which are applied at its ends 
in planes normal to its axis. 


6) 
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3*2 Shear S'taress s 

Since the shear stresses are developed by the 
torsional loading and radial shear forces (transverse 
shear), the shear stresses developed due to torsional 
loading are analysed as follows: 

■When the torque is applied to the shaft it 
twists, each section rotating about the longitudinal 
axis. Plane sections remain plane and radii remain 
straight. Thus at any point the shear stress "XrORQ 
developed on the plane of the section, the magnitude 
of this stress is proportional to the distance from the 
center of the section and its direction is perpendicular 
to the radius drawn through the point. Along with this 
shear stress there is an equal longitudinal shear stress 
on a radial plane, and equal tensile and compressive 
stresses at 45°, In addition to these stresses, the 
elongation of the outer fibers due to twisting causes 
longitudinal tension over the outer half of the section 
area and longitudinal compression over the inner half, 
the total tension and total compression being equal. 

The maximum longitudinal tensile stress occurs at the 
surface, the equal maximum longitudinal compressive 
stress occurs at the center* For the conditions here 
assumed, these longitudinal stresses are so small as 
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to be negligible in practically all cases 


[9]. 


Thus shear stress 'TtorQ ^ shaft due to 
torque T varies uniformly from zero at the center to 
maximum at the outside. The maximum shear stress takes 
place at the outer surface i.e. when r^ = d/2 and is 
given as 


T r. 


r 


(3.1) 


max “ J 

Since most of the shafts have solid/ circular cross- 

sections/ the polar rroment of inertia is equal tO/ 

,4 


J ^ 


tr d 

15 " 


(3,2) 


Thus maximum torsional shear stress for solid/ circular 
cross-sections becomes 




max 


16T 

TTd^ 


(3.3) 


A shear stress is also created by the shear force or 
transverse loading. First/ a shaft is assumed to be 
modelled as lengthwise elements. When the shaft is 
loaded/ the elements within the shaft try to slide on 
each other# but since these elements are rigidly attached 
a shear stress is developed. 


For any cross-section/ the maximum transverse 
shear stress occurs at the neutral axis# for circular 
cross-section the maximum transverse shear stress is 


given aS/ 
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~ 3 A (3.4) 

where/ V = Shear force 

A = Cross— sectional area. 

Since the transverse shear and torsional shear stress 
are in the same plane/ they can be combined into one 
stress. Assuming constant torque for a circular cross- 
section, the maximum torsional shear stress is at any 
location on tlie circmimf erence. However, the transverse 
shear stress ranges from zero at the outer fibers to a 
maxim\im at the neutral axis. The location of the neutral 
axis depends on the direction of the shear force, as 
shown in Pig. (3.1). Since the transverse shear is a 
maximum at the neutral axis, the limiting cases for 
shear stress occur at the outer fibers, as indicated in 
Fig. (3.1). Depending on the directions of the loads, 
the shear forces add on one side and subtract on the 
other. Since the shaft rotates ^each point on the outer 
surface experiences a fluctuating nonzero mean shear 
stress. 

3. 3 Normal Stress ; 

The normal stress is sometimes defined as the 
force per unit area. For a shaft, a normal stress is 
produced by an axial force and a bending moment. The 
nojrmal stress due to an axial force^^P , for a solid. 





b) Element A, shear stresses add. 

^'^Tr 


4 


i 



c) Element B, shear stresses subtract. 
Fig, 3.1 Location of maximum shear stress. 
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circular cross-section iS/ 


'Tm =■ 


(3.5) 


where/ 

= normal stress for axial load 
d s: diameter of cross-section 
D^hen normal stress due to bending is considered one 
side of the shaft is in tension, while the other side 
is in compression. The surface where stress is zero, 
is known as neutral axis* The maximum normal stress 
for bending for solid circular cross-section is 


<rmax 


32 M 
m d^ 


(3.6) 


Combining the two normal stresses is similar to the 


case with two shear stresses. If the axial force for 
the shaft is assumed not to vary, the forces produce 
a constant normal stress across the cross-section. At 
any cross-section, the normal stress developed by the 
bending moment is zero at the neutral axis to a maximum 
at the outer edges. For each cycle that the shaft 
rotates, any one point on the outer surface experiences 
a maximum compressive stress and a maximum tensile 
stress from the bending moment. Theresfore, the shaft 
has a non-zero mean fluctuating load. If the shaft 
is not axially loaded, the shaft develops a zero mean 
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fluctuating load [12]. The limitiiig cases are the two 
locations presented in Pig. ( 3 . 2 ). 

- • 4 Stress Concentrations s 

Failure generally is initiated at locations of 
geometrical discontinuities, which are known as stress 
concentrations. Ihese geometrical discontinuities, such 
as girooves, keyways, fillets, and abtupt changes, cause 
irregularities in the stress patterns. A stress concen- 
tration factor is defined as 

K = Actual stress at cross-section 
- Nominal stress at cross-section 

... (3.7) 

Many methods have been used to find these stress conc^i- 
tration factors for different types of geometry and 
loading. These values of stress concentration factors 
can be found from various references. For example. 

Fig, (3.3) to (3.5) are taken from Peterson [is], and 
mathematical equations are given by Roark [9]. 

Considerations must be given to stress concen- 
trations, since it is almost inpossible to design a 
shaft without them. Minimizing tiie effects of stress 
concentrations becomes a necessity for an acceptable 
design. Some ways to decrease stress concentration 
effects for shafts are to use large fillet radii, to 
add grooves and to place stress concentrations in low 
stressed areas. 



4-B 



a) Shaft with axial load, P and applied load,W. 




b) Element A, normal stresses add. 



c) Element B, normal stresses subtract. 
Fig. 3.2 Location of maximum normal stress. 






Fig. 3.3 Stress concentration factor K|^ for bending' of 
shaft with shoulder fillet. Ref. [.13J. 
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Fig 3.4 Stress concentration factor for torsion of 
shaft with shoulder fillet. Ref.[l3j 
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Effects of the stress concentrations can be 
quickly evaluated with a computer program* Also, a 
major advantage of a program is its ability to show 
locations of low stress for placement of stress concen ' 
trations* 


3*5 Failure Theories? 


Before introducing the failure theories, an 
expression must be developed to combine a steady and 
alternating stress to an equivalent steady stress. There 
have been many equations developed for test data to 
relate non-zero mean alternating stress to an equivalent 
zero mean alternating stress. The most conservative 
equation to perform this transformation is a Soderberg 
equation [3,lo]. 


For normal stress, the Soderberg equation is 
O' = (Tav '3-8) 


where. 


av 
tTr 
yP 
OT e 


= average normal '.stress 
= 1/2 normal stress range 
= materials yield stress 

= materials endxirance limit for zero mean 
stress loading 


K « stress concentration factor 
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For shear stresses, the Soderberg equation becomes 

-r . X + Xr (3.9) 

(Te 

where shear stress 

"^j,= 1/2 shear stress range 
K^=: stress concentration factor for shear 

Equations (3.8) and (3.9) take into account stress con- 
centrations. However, the stress concentration factors 
are usually applied only to the alternating stresses. 

The stress concentrations in the two equations are 
different. For example in equation (3.3), K is for 
bending, whereas in equation (3.©) is for a torsiona] 
load. Since a rotating shaft is under consideration, 
the average normal stress for equation (3.8) is equal 
to the stress due to the axial loading. The amplitude 
of the alternating nomal stress is picduced by the 
bending moment and by any variations in the axial load- 
ing. However, the bendiing stress does not have to be 
completely reversed, but for the worst possible design 
case, the bending stress would be totally reversible. 

For equation (3.9), the average shear stress is due to 
the constant torque. The amplitude of the alternating 
shear stress is equal to the transverse shear or to any 
fluctuation in the torsional shear stress. Again, for 
the worst possible case, the transverse shear is presumed 
to be completely reversed. 



Since most shafts are made of ductile material o,- 
the two failure theories most accurately representing 
them are the maximum shear stress and distorsion energy’ 
failure theories. The distorsion energy failure theory 
fits experimental results for ductile materials extremely 
well [ll] and [s], 

For stresses in two dimensions, the theory gives 


g2 



^•gT y 





xy 


fs.io) 


where S = equivalent normal stress 

^ = normal stress in x— direction 

V yc 

s normal stress in y--dir action 

rf" « shear stress in X/y plane# 
cxy 

For the case of rotating shafts equation (3# 10) rednr^ 


to 

s2 ^^2 ^3-^2 (3.11) 

After substituting equations (3.S) and (3.9) into the 
equation (3.11) and taking the square root, equation 


(3.11) becomes. 




e 

(3.12) 


The working stress for uniaxial loading can be defined 


as 


(TyP / FS 


(3-13) 



where 
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0^ = Working stress 
FS= Factor of safety 

After substituting the working stress/ (T * fot the equi 
Valent stress, 3/ equation (3.13) is equal to 



The maximum shear stress failure theory is also used 
for predicting failure for ductile materials/ although 
it is less accurate than the distorsion energy theory. 
The first step is to find the maximum shear stress for 
a given loading. For two plane loading/ the maximum 
shear stress is 



For shafting, this reduces to 


(3.15) 


v' (f )2 + Y" 


(3.16) 


After substituting equations (3.8) and (3.9) into 
equation (3,16), the maximum shear stress theory predicts 


failure as. 



(3.17) 
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For pure shear loading, the working stress is 

max ~ *T*yP / FS (3. 18) 

where, 

= Shear yield stress 
FS = Factor of safety 

The shear yield stress for uniaxial loading is 


'TyP ~ ^P/2 


(3.19) 


By substituting equation (3.19) into equation (3.18), f 


working stress becomes 


T 

\ tr 


(3.20) 


Finally, after substituting the equation (3.20) into 
equation (3.17), equation (3.20) becomes 


0.5 


= ,/ 1(0- +^r)2 + (r 


4' av ^ 




(3.21) 

Depending on which type of failure theory is to be used 
the factor of safety can be determined from equation 
(3.14) or equation (3.21). For both failure theories, 
it is assumed the normal and shear stress maximums 
occur simultaneously, which is not necessarily true. 
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3.6 

Madiine members are often found to have failed 
under the action of repeated or fluctuating stresses# 
though the analysis reveals that the actual maximum 
stresses were below the ultimate strength of the material 
and below yield strength. This type of failure occur 
because stresses are repeated a very large number of 
times and this failure is called a fatigue failure. 
Fatigue failures mainly occur at change in cixiss- section# 
keyv-ray^or a hole. 

When a fatigue strength is to be considered 
in design and analysis of rotating shaft the modified 
endurance limit of that shaft is calculated [lo], as 
given in equation. 

S = K K. K S' (3.22) 

e a b c e 

S = modified endtarance limit of mec±ianical 
e 

element 

s' = endurance limit of rotating beam 
e 

specimen 

K = surface factor 

a 

= size factor 

K IS modifying factor for stress concentration 

c 



where# 



58 


By itiodified endurance limit of a material is usually 
meant that the maximiim stress which can be reversed on 
large number of times without producing fracture. For 
steels the endurance limit of rotating beam specimen is 
calculated by equation, 

s; = 0.5 (3.23) 

where, 

^ut = ultimate tensile strength. 

The factors , Kj^ and are governed by equation 
surface finish factor is given in [lo] as follows! 

K =0.73 for ground starface finish 

(3.24) 

“ 0.84 for polished surface finish' 

Size factor is taken into account since when alternating 
stresses are present in bending and torsion the diameter 
of shaft has also some effect on fatigue failure [2,10], 
it is given as, 

= 1.189 (3.25) 


where, 

d = diameter of shaft in mm 
modified factor for stress concentration is given as 


where 


K 


c 



(3.26) 



fatigue stress concentration factor 
which is given by equation (3.27) 
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= 1 + q (K^ - 1) (3.27) 

where, 

q = notch sensitivitif 

= theoretical stress concentration factor 

notch sensitivity q [l3,16] can be found from Pig. (3.6), 


When this modified endurance limit is calcu- 
lated from eq.(3.22)/ this modified endurance limit 
is used in failure theories discussed in Section 3*5 
when fatigue strength effect is to be considered. 



NOTCH SENStTIVITY 



Fig 3.6 Notch sensitivity, q Ref. [16], 



CHAPTER - 4 


DESIGN AND ANALYSIS OP ^^ULTI~SPEED GEAR BOX SHAFTS 
4* 1 Introduction : 

This cihapter gives a brief idea about/ how to 
analyse gear box shafts. As discussed in Chapters 2 
and 3/ to analyse -any shaft# loading on the shaft 
should be known. Since gear box provide for a wide 
range of cutting speeds and torques from a constant 
speed power input to get proper speeds and torques at 
the spindle*. The loading on these shafts are mainly 
torques# tangential and ’radial point loads when the 
gears on these shafts mesh wiih. gears on other shafts. 
Bending moments and shear forces on these shafts due to 
the tangential and radial point loads can be calculated. 
Then these torques# bendir^ moments and shear forces on 
the shafts are taken for the analysis of shafts. As 
discussed in Chapter 3 # using these torques# bending 

moments and shear forces# factor of safety at any point 
along the shaft is calculated. Maximum deflection of 
the shaft is calculated since this is also a factor used 
to analyse a shaft. 
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■^•2 Force Analysis; 

Since iTailt.i~speed gear box gives a wide range 
of speeds and torques at spindle, the main loadings 
on the intermediate shafts are torques, tangential and 
radial loads. In this present work the gear box consis- 
ting of only spur gears (which are most common), having 
same module are considered. However, since there is 
a provision to consider axial loading in the program 
developed, the gear box shafts having worm or other 
gears can also be analysed by adding a new subroutine 
for their force analysis. When two gears mesh 
each other, forces are exerted on each gear along 
the pressure angle 0, where the two pitch circles of 
these o gears meet as shown in Fig. 4.1(a). Consider- 
ing free body diagram Fig. 4.1(b) it can be seen that 
the force on the pitch circle of gear can be replaced 
by a torque and a point load called normal force on the 
shaft. Thus proceeding this way we can find the torque 
and point loads on the shaft at different places where 
the gears on this shaft meshes with the gears on other 
shafts. 


F 

P 


n' 

21 

n 

12 


In Fig. (4.1) the forces are defined as 

normal force exerted by gear 2 against gear 1 

normal force exerted by gear 1 against gear ? 
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^21 ~ tangential force on shaft (in horizontal plane) 

^21 ~ ^a<3ial force on shaft (in vertical plane) 

- torque on gear 2. 

When these tangential and radial forces on the shaft 
are calculated # the bending mori->ents and shear forces 
in horizontal and vertical planes are calculated, then 
these are combined together and considering theories 
of failure as discussed in Chapter 3 > the factor of 
safety is calculated. It is seen that, when the shafts 
iagear box are not coplaner, then also the values of 
resultant bending moments and shear forces on the shafts- 
donot change. This is due to the fact that, the tangen- 
tial and radial forces in two planes found at any loca- 
tion are governed by 


Similarly the resultant bending moments and shear forces 
are governed by 


(bm)^' + (BM)^ 

and - 

(SP)^ + '(SP)^ 
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where, 

- resultant bending moment. 

(BM)pj - bending moment in horizontal plane. 

(BM)^ ~ bending moment in vertical plane. 

(SF)r - resultant shear force. 

(sF)r. - shear force in horizontal plane. 

(SF)^ - shear force in vertical plane. 

Thus for a compact gear box the shafts can be noncoplaner. 
This can be seen from Pig. (4.2). In the program deve- 
loped the equations used to calculate the forces in 
horizontal and vertical planes from normal force on the 
gears are as 

F„ = cos (110° - e) (4.4) 

n 

F^ = f’^ sin (ilO° - 0) (4.5) 

where, 

9 — a3ngle between two shafts with respect 

to horizontal plane as shown in Fig, 

(4.2 ) (c) - 

The angle 110° in equations (4.4) and (4.5) is the 
addition of pressure angle and the angl^ botveen the 
tangent to pitch circles and line joining the center 

shown in Fig. (4.2) 


of the gears as 
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^ ?2gsiqn and Analysis of Gear Bo.; .Shafts t 

Before designing the gear loox shafts the input 
horse pov;er and speed of the motor should be known. By 
using the method as described in [is] the gear box is 
designed, from which the inputs to design and analyse 
the gear box shafts taken are number of shafts, number 
of gears on each shaft, number of teeth on each gear, 
module, pressure angle and speed ratio. In the present 
developed program it was not possible to get automati- 
cally the different possible arrangement of the gears 
and the sequence of the gears in which they mesh to get 
possible different output speeds in a particular arrange- 
ment, Thus these sequence of the gears which mesh to 
give proper output speeds is also one of the inputs and 
this can be done for number of possible different arrange' 
mont of the gears. Some rules are to be followed to get 
a proper design and analysis of gear box shafts, and 
also to give proper inputs to the program such as the 
different possible gear arrangements and the placement 
of the gears in a particular gear arrangement. These 
rules are as follows, [ 3 ] 

a) One set of gears must be completely disengaged 
before the other set begins to mesh. 

The sum of the teeth of mating gears in a given 
stage must be the same for same module in a 


b) 



6 S 

clustsred set since the i_c;nter distance must 
remain same. 

c) The minimum .ixid grp between Wo e^jaccnt gear':; 
must be at l-^a'c tvrni gc'"'r x.'dths. 

d) The clearance between the fixed gear and 
sliding gear in the small and medium machine 



toolv gear box should be 1*5 ix> 3 rnm# 


The diameter of 

a gear box shaft is calculated 

roughly 

by usir^ 

eg. (4 

.6)^ [ 3 ] 



^(I) 


(4.6) 

where 





'^(l) 

^(I) 

*thi 

= diameter of I shaft 

= maximu. torque on the I shaft 



= permissible shear stress. 



This rough diameter of a shaft is standardised 
from a table 4.1 [l6]. 

Then if the hey\-;:ys are present their locations 
and the steps on the shaf .:s on which the gears are 
mounted can be given as input if necessary. After giv- 
ing the bearing positions, the shaft analysis starts. 
For the shaft analysis as discussed in Section (4.2) 
the forces and torques are calculated then these loads 
are used to calculate bending noments and shear forces 



TABLE 4.1 


Table to select the standard diameter of the shaft 


Standard 
Diameter 
in cms. 

Standard ; 
Diameter 
in cms. 

standard ! 
Diameter • 
in cms. j 

Standard 5 
Diameter j 
in cms . | 

standard 

Diameter 
in cms. 

0.6 

0.7 

0.3 

T 

{ 

1 

0.9 

1.0 

1.1 

1.2 

1.4 

1.5 

1#6 

1.7 

1.8 

2 .0 

2.2 

2.4 

2.5 

2.6 

2.3 

3.0 

3.2 

3.5 

3.6 1 

3.3 

1 

1 4.0 

4.2 

4.5 

4.8 

5.0 

5.2 

5.5 

5.6 

5.8 

6.0 

6.2 

6.5 

6.3 

i 7.0 

i 

7.2 

7.5 

7.3 

! 

3.0 

i 

1 8.5 

! 

X 

9.0 

9 *5 

1 10.0 

j 


Using any failure theories as discussed in Qiapter 3 . 

the factor of safety and the maximum deflection at any 
point for a particular gear sequence which has maximum 
torque and bending moment in particular gear arrange- 
ment are calculated. Thus a shaft having least lengtii 
and safest one is chosen as the best shaft out of the 
possible options given. 



CHAPTER-5 


PROGRAMMING CONSIDERATIONS AND EXAI^LES 

5.1 Computer Aided In-beracbive Design and Analysis ; 

Computer aided design is a technique in which 
man and machine are blended into a problem ~ solving 
team, coupling the best characteristics of each. The 
result of this combination works better than either man 
or machine would work alone, [l?] 

Computer aided design and analysis consists of 
five stages as, design logic, equations and computations, 
design checking, engineering paper work generation and 
graphical layout of design. 

In design logic the computer's decision making 
capability can be utilized in selecting the best of a 
niimber of possible design alternatives. In the present 
work the program can select the best design after analys- 
ing each alternative designs. 

Almost all the equations and computations can 
be put in d form s uitable for solution. In the present 
work the transfer matrix method is programmed for itera- 
tive calculations, for analysis of shafts as discussed 
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in chapter 2 . Design checking can be a time constymlhg 
manual task, especially if n'umerous critical calculations 
are involved. Thus the program developed here takes care 
of this basic function, with appropriate data substitu- 
tions for specific design checks, and thus after analys- 
ing it selects the best of all. 

The computer is ideally suited for use in gene- 
rating much of the engineering paper work necessary to 
the successful implementation of a design. 

Graphical layout of the designed product can be 
seen on graphics terminals* The graphical plots of 
designed and analysed product can also be seen and from 
which the designer can decide vlaether the design is safe 

or not. 

The term interactive design and analysis means, 
the designer has some control over the intermediate results 
as well as final results, while the program is in execu- 
tion mode. Designer is provided with the facility of 
selecting the required options, and/or sometimes allowed 
to chairre the values of parameters calculated. All this 
logic is so developed that the program will automatically 
tahe care of the changes to he made according to the new 
information given by user. Ihus interactiveness of the 
program is a very powerful tool. 
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'The present program has been v/ritten in Fortran-^ 10, 
developed and tested on DFC-1090 system. In the present 
chapter the major programming considerations are put in 
the form of flow charts, also the purpose of different 
subroutines is explained, 

5,2 The Main Programs 

The purpose of main program is to accept the 
input data interactively to calculate some of the para- 
meters in itself and others by calling subroutines and 
lastly to print the results. Input data is stored in 
one file, which contains answers for the interactive 
questions which are variables, or single dimensional 
arrays. The output of -the program is stored is one 
separate file which can be printed to get a permanent 
record and the same file can be used as input to tiae 
graphics program. 

The program is written in a very general way 
so that it can handle number of different cases. Most 
of the interaction dialogues are built in the main 
program. This program developed can analyse shafts 
having any type of loadings and circular cross-section. 

It can also design and analyse the gear box shafts cal- 
culating the forces and torques on it, where input para- 
meters are only its constructional, gears and input 

power# 
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VJliile developing the pirogram it is kept in mind 
that the memoary locations are as less as possible. Some 
3 D matrices are used in the program these are very 
limited and the 3rd dimension o£ these matrices is to 
lable# just to remember these matrices. 

Separate subroutines are written for purposes 
described here. Plow cdaarts for some subroutines are 
given. 

1 , Subroutine MATMUt 

Purpose - To multiply the 2-D matrices 

2 , Subroutine AMAT 

Purpose " To create an identity matrix 

3, Subroutine APINS 

Purpose — To create an in— span point matrix 

4, Subroutine - MAT 3D 

Purpose - To multiply the 3-D matrices 

5 , Subroutine MAT 2 3D 

Purpose - TO multiply the 3-D and 2-D matrices 

6 , subioutine FORC 

Purpose - To calculate the forces on shafts 
in a gfear boX/ when given gear box 
constructional parameters and input 
power ana sequence of gears as ir^ut 
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7 . Subroutine ANLYS 

Purpose - To analyse the shaft ^calculation of 
stress concentrations at any change 
in cross-section and at keyways# 
consideration of fatigue failure is 
included in this subroutine* Finally 
using all these values and forces on 
the shaft considering distorsion 
energy or maximum shear stress failure 
theories, normal stresses and factor 
of safeties are calculated along th.e 
shaft. 

5.3 Graphics Program t 

The graphics program is developed by usirg 

General Purpose Graphics System manual L2o] and [iP]. 

of 

Once the main program is executed the output^/this program 
is used as input to the graphics program. The purpose 
of this program is to draw the deflection, slope, bend- 
ing moment, shear force, stress distribution and factor 
of safety plots along the shaft. Since from these plots 
the designer will have quick idea about how safe the 

shaft is. 

The interactiveness of the program can help 

select any plot ■which he wants to see. 


the designer to 
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There Is also a provision to draw all these plots on the 
same screen. ^Ihen shaft selection for particular purpose 
is to he made then these plots help the designer to quite 
an extent. The problems like giving window for each 
plot as input, etc. are overcome in the present developed 
graphics program. Since this program automatically scales 
the screen accordingly oncce by reading the input, by 
calculating the maximum and minimum of the values. 

The flow chart of tie graphics program is shown 


here. 



MAIN PROGRAM 


START 


Call 

subroutine 


/ gear box^ 
^ shaft is to 
'^xanalysed ^ 


Input geometric propearties 
material properties and 
boundary conditions 


Input loadings# locations 
of the loading in vertical 
and horizontal planes 


10 

#2 

(planes) 


Determine no . of 
Ideations of 
changes 


Locate positions 
of transfer matrices 


,, . - i 

Initialise subglobal 
matrix to identify 
matrix call subroutine 
AMAT 


Create field matrix 
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( END } 





SUBROUTINE PORC 


Prom Main Program 


Input HP , no , of g ears , no . of teeth 
on each gear, module, pressure angle 
and location of these gears on each . 
■ shaft / 


Find riax. 
torque on 
each 
shafts 


/^If \ 
shaft is 
to be 
(designed 


Find diameters 
of shafts 
using torque 
equation 


Standardise 
the diameter 
from the 
table given 
on the 
terminal 


Give staps 
on the shaft 

if necessary 
giving filled 
radii 


Find the forces on 
each shaft other quan- 
tities whicii will be 
used by the main 
program directly 


return 


To main program 







SUBROUTINE ANLYS 


From Main Program 


Input the yield stress , ultimate 
tensile strength and endurance 
limit for the shaft material i 


( Input 
stress 
concen- 
tration; 


to give stress 
concentrations for stuped 
'-v. portion and keyways?,---^ 


Input fillet radius at 
stepped portion and 
length of keyway s etc. 


Calculate stress cocentrations 


Select any one failure theory 
I among dis torsion energy and 
maximum shear stress theories 


Find factor of safety and 
I normal stress 


Output the results 


Output smallest factor of safety 
and largest deflection along the 
shaft — i 


RETURN 


To Main Program 











GRAPHICS PROGRAM 


START 


Input values of deflection, 
slope bending moment, shear 
force stress distribution, 
factor of safety 


Input about shaft 


Calculate the maximum and minimum 
of the input values of which plots 
^re to be drawn j 


Scale the windows according to their 
max, and min. values to get plots on 
the screen 


Interactive features to dravr any 
plot on the screen 


Plot on the screen by 
using (RPGS routines 
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^ ^ .'gy ample 1 : 

Tliis is a analysis problem. The aim of tihis 
problem is to got normal stresses and factor of safeties 
along the shaft. The geometric p 3 Xipertios and the load- 
ings are as shown in Pig. (5.1)/ all the unit in this 
problem are British units. 

The inputs for this example at the time of the 
execution of the program are given on the terminal as/ 

DO YOU WANT TO DESIGN AND ANALYSE GEAR BOX SHAFTS? TYPE 
Y/N : N 

DO YOU WAI^IT TO DESIGN A SHAFT? TYPE Y/N s N 

PLEASE TYPE IN THE VALUES FOR YOUNG'S MODULUS AND SHEAR 

MODULUS; 30000000.0, 11000000.0 

For boundary conditions lihe free end type 0, for pinned 
end type 1, and for fixed end type 2. 

PLEASE TYPE IN THE LEFT AND RIGHT END BOUNDARY CONDITIONS 
and NUMBER OF SUBGLOBAL MATRICES TO BE F0RI'4ED : 1, 0, 2 
Type in the geometric properties and loadings for 
Horizontal plane. 

PLEASE TYPE IN NUMBER OF FIELD MATRICES IN EACH SUBGLOBAL 
MATRIX TO BE CREATED : 7/ 2 

PLEASE TYPE IN NUMBER OF CONCENTRATED LOAD OCCURRANCES 
IN EACH SUBGLOBAL MATRIX ; 4, 0 

IS THERE ANY AXIAL FORCE PRESENT? TYPE Y/N ; N 
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PLEA3E TYPE IN DIAMETERS OF SHAT’T IN EACH FIELD MATRIX? 

0 . 6 , 0 . 6 , 0 . 6 , 0 . 6 , 0 . 8 , 0 . 3 , 0 . 3 , 0 . 3 , 0.8 

PLEASE TYPE IN SPAN OP EACH FIELD MATRIX ? 2.5, 0.5, 

2.0, 1.0, 2.0, ?,0, 4.0, 1.0, 1.0 

PLEASE TYPE IN MlAGNITUDES OF DISTRIBUTED LOADS IN EACH 
FIELD MATRIX : 0.0, 0,0, 0,0, 100.0, 100.0, 100.0, 100,0, 

100 . 0 , 0.0 

PLEASE TYPE IN MAGNITUDES OP GRADIANT OF LINEARLY 
DISTRIBUTED LOAD IN EACH FIELD MATRIX ? 0.0 0.0 0,0 

0.0 0,0 0.0 0,0 0.0 0.0 

PLEASE TYPE IN DISTANCES OF CONCENTRATED LOAD OCCURRANCE 
FROM LEFT END : 2,5, 3.0, 3.0, 11.0 

PLEASE TYPE IN THE MAGNITUDES OP SHE-AR FORCES ? -200.0, 

- 200 . 0 . 

PLEASE TYPE IN THE MAGNITUDES OP CDNC®7TRATED MOMENTS? 200 .0 
PLEASE TYPE IN THE MAGNITUDES OP SPRING CONSTANTS OF 
TRANSLATIONAL SPRINGS ? 5000000.0 

In simi-lan "typs in tLs data fon Vsirtical plane. 

Input for the Torque are, 

Giv^ the looundary conditions, if shaft is fixed as 1, 
if free or pinned as 0. 

PLEASE TYPE IN THE BOUNDARY CONDITION AT LEFT AND 
RIGHT END ? 0, 0 

PLEASE TYPE IN THE NUI4BER OF TORQUES APPLIED ? 3 
PLEASE TYPE IN THE MAGNITUDE OP EACH TORQUE ? 600.0, 

600.0 , -1200.0 
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PL3ASE TYPE IN THE LOCATIONS OF TORQUES FRONi LEFT 
END s 0.0, 10.0, 15,0 

DO YOU WANT TO USE DISTORSION ENERGY FILURE TNEORY 
TYPE Y/N j Y 

PLEASE TYPE IN THE YIELD STRESS AND SI'IDURANCE LIMIT 

OP SHAFT MATERIAL ; 130000,0, 120000.0 

DO YOU WANT TO GIVE STRESS CONCENTRATION FACTORS FOR 

BENDING AN TORQUE? TYPE Y/N ; N 

PLEASE TYPE IN THE FILLET RADIUS ; 0.01 

DO YOU WANT TO CONSIDER FATIGUE FAILURE AT THIS FILLET? 
TYPE Y/N : Y 

PLEASE TYPE IN THE NOTCH SENSITIVITY s 0.9 

PLEASE TYPE IN NUMBER OF KEYWAY S PRESENT, IP NO KEY-WAYS 

TYPE 0 ; 0 

The results of the analysis are stored in DATA files, 
and are given in the following tables (5.1) to (5,4). 
These tables are used by the graphics program and the 
plots of the different values can be seen on the graphics 
terminal . 

smallest . 

The/factor of safety is 2.3676 at 2.39 inches 

from left end. The maximum deflection is 0.02916 inches 

at 3.57 inches from left end. 
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b) t-oading in vertical plane. 


Fig 5.1 Loadings for Example 1 . 
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5.5 Example 2; 

is a design and analysis problem. The 
aim of this 'oroblem is to get safe diameter of the shaft 
for a given loading, considering stress concentrations 
at keyways and to find normal stresses and factors of 
safety/. The geometric properties and the loadings are 
as shovjn in Pig. (5.2). The problem is solved in SI units. 
The input is as follows, 

DO YOU WANT TO DESIGN AND ANALYSE GEAR BOX SHAFTS? 

TYPE Y/N : N 

DO YOU WANT TO DESIGN THE SHAFT? TYPE Y/N : Y 
PLEASE type in THE MAXIMUM TORQUE, BSNDItIG MOMENT ANL 
PERMISSIBLE SHEAR STRESS s 4792.0, 4500.0, 420.0 
The calculated diameter of the shaft is 3.8 cm, 

DO YOU WANT TO STANDARDISED THE DIAMETER OP THE SliAFT? 

TYPE Y/N : Y 

The table of standard diameters of shaft (as shown in 
chapter 4) is shown on the Terminal. 

PLEASE TYPE IN THE NEAREST STANDARD DIAMETER IN QiSsi.O 
PLEASE TYPE IN THE VALUES FOR YOUNG'S MODULUS AND 
SHEAR MODULUS s 21000000.0 , 8400000.0 

Type in the locations at the loadings and their magnitudes. 
PLEASE TYPE IN THE LEFT AND RIGHT END BOUNDARY ODNDITION 
AND NUMBER SUBGLOBAL MATRICES TO BE FORMED : 0, O, 3 
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PLEASE TYPE IM THE ^7U1'IBER OP FIELD MATRICES IN EACH 
SUBGLOBAL MATRIX TO BE CREATED s 1 , 1 , 1 

PLEASE TYPE IN THE NUMBEP OP CONCENTRATED LOAD OCOJRRANCES 

IN EACH SUBGLOBAL MATRIX r 1, 0, 1 

IS THERE ANY AlCIAL FORCE PRESENT? TYPE Y/N ; N 

PLEASE TYPE IN SPANS OP EACH FIELD MATRIX ; 15.0, SO.O, 

25.0 

PLEASE TYPE IN THE DISTANCES OP CONCENTRATED LOAD 
OCaJRRANCE FROM LEFT END ; 0.0, 120.0 

PLEASE TYPE IN THE MAGNITUDES OF SHEAR FORCES s 300,5, 
160.3 

Inputs for Torque - 

PLEASE TYPE IN THE BOUNDARY CONDITIONS AT LEFT AND 
RIGHT END : 0, 0 

PLEASE TYPE IN THE NUMBER OF TORQUES APPLIED t 2 
PLEASE TYPE IN THE MAGNITUDES OF EJvCH TORQUE ; 1792.0 

-1792.0 

PLEASE TYPE IN THE LOCATIONS OF TORQUES PROM LEFT END. 

0 . 0 , 120.0 

DO YOU WANT TO USE DISTORSION ENERGY FILURE THEORY? 

TYPE Y/N : Y 

PLEASE TYPE IN THE YIELD STRESS AND ENDURANCE LIMIT; 
3164.0, 2953.0 

PLEASE TYPE IN THE NUMBER OF KEY-VJAYS PRESENT ; 2 




Fig,5.2 Example 2 . 




Distance along the shaft cm ^ 

Fiq 5.3 Riots for shaft of Example 2. 
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PLEASE TYPE IN ItiE STRESS CONCENTRATION AT EACH 
KEY vjaY : 1.5, 1.5 

PLEASE TYPE IN TtiE LEFT STARTING POSITION OF E/^CH 
KEY WAY s 0.0, 110.0 

PLEASE TYPE IN THE LENGTH OP EACH KEY W7Y s 10.0, 10.0 
The results of this designed and analysed shaft are stored 
in the DATTi. FILES, ihe plots of the resultant values 
nomal stresses and facto res of safety are as shown in 
the Pig. (5,3). 

The diameter of the designed shaft is 4,0 cm. 

The smallest factor of safety is 3.92 which is at 15.0 cm 
from left end. The maximum deflection is 0.19 cm at 
120.0 cm from left end i.e, at right end. 

5 . 6 Example 3 s 

This example illustrates a 4~speed machine tool 
gear box. Ihc problem aims to design and analyse the 
gear lyox sha-^^ts for a given gear arrangement. This 
problem is executed for different possible gear arrange- 
mc.^nts in the gear ^oox. Then a shortest and an appropriate 
gear box shaft is selected from the different specified 
options. Tw gear arrangements of a 4-speed gear box 
are shown in Fig. (5.4)a and (5.4)b. Tho shortest and 
appropriate gear box shafts are for gear arrangement 2 
as in Fig. (5.4)b. Bie Inputs for the gear arrangement 2 
(taKen from [ 3 ]) are as folXows=(all Inputs are in SI Units) 
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DO YOU want to design AND ANALYSE GEAR BOX SHAFTS? 
yes Y/N ; Y 

PLE-AjE type in the module of the gears in CMS s 0.3 cm 
PLEASE TYPE IN THE PRESSURE ANGLE OF THE GEARS IN 

degrees : 20° 

PLEASE TYPE IN THE HORSE POWER OF THE INPUT MOTOR; 10.0 
PLEASE TYPE IN THE SPEED OP THE INPUT MOTOR; 

PLEASE TYPE IN THE NUMBER OP SHAFTS IN THE GEAR BOX; 3 
Give code 1 for the sliding shafts and 0 for non-sliding 
shafts. 

PLEASE TYPE IN THE CODES FOR SLIDING AND NON-SLIDING 
SHAFTS ; 1, 0, 1 

PJjEASE TYPE IN THE ANGLE AT WHICH THE SECOND TO LAST 
.SHAFTS A,RE# WITH THE FIRST SHAFT ; 0.0, 0.0, 0.0 
PLEASE type in THE DENSITY, BREADTH OP GEARS : 0.0072,3.0, 
PI, EASE TYPE IN YOUNG'S MODULUS AND SHEAR MODULUS FOR 
•SHAFT MATERIAL ; 210000000.0 , 8400000.0 
PLEASE TYPE IN THE NUMBER OP GEARS ON EACH SHAFT; 2, 4, 2 

PLEASE TYPE IN THE NUMBER OP TEETEI ON SACK GEAR ON EACH 

SHAFT : 24, 20 

PLEASE TYPE IN THE POSITIONS OP THE GEARS FROM LEFT 
end BEARING ON EACH SHAFT; 2. '0, 5.0, 14,0, 17,0 
PLEASE TYPE IN THE POSITION OF THE RIGHT END BE7.RING ; 'i9,o 

PLEASE TYPE IN THE SEQUENCE OP TIE GEARS FOR PARTICULAR 

COMBINATION ; 1 ,'l ) 2 , 1 
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The inputs tor gear arrangement 1, are also given in the 
similar manner. 

The diameters of all the shafts are calculated 
and then they are standardised from the table typed 
on the terriiinal of the time of execution. The steps 
on the shafts and the lengths of the keyways are also 
given. Then the analysis of the intermediate non-sliding 
shafts start. 

From the plots of these results of 'the two 
arrangements shown in Fig, (5.5)a and (5.5)b it can 
be seen that^for arrangemest 2 tlae intermediate non- 
sliding shaft design is safe and -the length of the gear 
box is small. 

The diameters of the shafts in gear arrangements 1 
and 2 are same, as same maximum torques on each shafts 
are present, only the lengths of the shafts in gear 
arrangGraent2are smaller than the gear arrangement 1. 

The diameters of the shafts obtained from the 

design, in gear arrangement 2 are, 

diameter of shaft 1 is 2.0 cm. 

diameter of shaft 2 is 2.5 cm. 

diameter of shaft 3 is 3.0 cm. 

Lengths of the shafts in gear arrangement 1 

are 25 cm. 

■ Lengths of the shafts ^ arrangement 2 

are 19 cm. 






Fig 5.4 (a) Examples. 4 Speed gear-box 1^^ arrangement. 







. J'Cl 

Fig 5.4(b) Example 3 , 4 Speed gear-box 2 arrangement. 





Resultant Shear force 



Fig 5.5(b) Plots for Examples 2 gear arran^ment 
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For shaft 2 in gear arrangement 1 the values of minimum 
factor of safety and maximum deflection are 1.749 and 
0.0016 cm, respectively. 

For shal't 2 in gear arrangement 1 tlie values of minimum, 
factor of safety and maximum deflection are 2.'''»25 and 
0.00075 cm respectively, 

5.7 E>:ample 4'; 


The aim of this problem is to design and analyse 
the IS speed gear box shafts/ for a given arrangement 
of gears. For a gear arrangement shown in Fig, (5,6) 
short and safe ^afts are designed and analysed. 

The inputs for this gear arrangement (taken from [ 15/18]) 

arc as follows: 

BO YOU WANT TO DESIGN AND ANALYSE GEAR BOX SH-FTS? 

TYPE Y/N ; Y 

PLEASE TYPE IN THE MODULE OF THE GEARS IN QIS : 0.3 on. 
PLEASE TYPE IN THE PRESSURE ANGLE OF THE GEARS IN 


DEGRESS : 20 

PLEASE TYPE IN THE HORSE POWER OF THE INPUT MOTOR: 6.711 
PLEASE WPS ™ 

please TWE in IHE cedes for eliding and HON-gLIDlNG 


SHAFTS : 1/ 0# ^ ^ ^ 

PLEASE TYPE IN THE ANGLE AT WHICH THE 
SHAFTS ARE WITH THE FIRST SHAFT : 0.0 


SECOND TO LAST 

, 0.0/ 0.0/ O.O/O.O 
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PLEASE TYPE IN THE VALUES FOR DENSITY AND BREADTH OF 
GEARS ; 0.0072, 3.0 

PLEASE TYPE IN THE VALUES FOR YOUNG'S MODULUS i\ND SHEAR 


MODULUS ; 21000000.0 8400000.0 

PLEASE TYPE IN THE NUMBER OF GE.ARS ON E7-GT SHAFT s 3, 6 , 4 , 3, 2 
PLEASE TYPE IN THE NUMBER OF TEETH ON EACH GUARS ON 
EACH SHAFT ; 20, 34, 47, 56, 43, 27, 31, 33, 34, 37, 

35, 32, 34, 20, 26, 35, 54, 47. 

PLEASE TYPE IN THE POSITIONS OF THE GEARS PROM LEFT END 
BEARING ON INTER14EDIATE SHAFT: 2.0, 11.0, 20.0, 2 3.0, 

32.0, 41.0 cms. 

PLEASE TYPE IN THE POSITIONS OF THE RIGHT END BEARING J 46 cm, 
PLEASE T£PE IN THE SEQUENCE OF THE GEARS FOR A PARTICULAR 
COMBINATION; 1, 1, 4, t, 4, 4, 1, 1 

The diameters of the shafts for this arrangement 


are calculated and 

PLEASE TYPE IN THE 
2.5, 2.6, 2.3, 3.0 
PLEASE TYPE IN THE 


typed 

ST7UTDARDISED DIAi4ETERS OF SHA.FTSs 
STEPPED DIAMETERS ON THE NON-SLIDING 


SHAiFTS: 2.5, 2.6, 2.8, 2.6, 2.5 

PLEASE TYPE IN THE NUMBER OF KEYWAY S ON EA.CH SHAFT; 5 
PLEA.SE TYPE IN THE POSITIONS OF THE KETNAYS ON THE 
SHAFT :0. 5 , 3.5, 13,5 , 30.5 , 42.5 

PLElSE TYPE IN IHE LENGTH OF E.’.CH KEfWAY = 3.0, 6.0, 
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PLEASE TYPE IN THE STRESS CONCENTRATIONS AT KEYWAYSsl.2 
PLEASE TYPE IN THE NOTCH SENSITIVITY FOR FATIGUE FAILUREsO.S 

The diameters of all the shafts are calculated 
and they are standardised from the table typed on the 
terminal, 

Ihe diameter of the shafts are as follows / 

The standardised diameter of shaft 1 is s 1.7 cm 

2 is s 2.5 cm 

3 is s 2.6 cm 

4 is s 2 .8 cm 

5 is s 3.0 cm 

Length of the gear box is 46 cm. 

The results of the analysed intermediate shafts are tabulated 
^DATA files. The plots for the 2na and -lao nonslldlhg 
shafts are shovin In Fig. (S.VVand Fig. (5.7 )b. 

Factor of safety for shaft 2 is 2.12 and for 

shaft 4 is 4.06 

Maximum deflection for shaft 2 is 0.0094 cm and 

for shaft 4 is 0*003 cm* 




Resultant ShcOir force 



•^0 • goo. ^00* loo. 

— UJD UOU33\|3(] 



Resultant Shear force 
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CONCLUSIONS 

6 . 1 Technical Summary : 

Computer aided design is a design process in 
which engineer uses the computer to design machine 
elements with the help of interactive computer programs. 
The present work is an attempt to develop an interactive 
design and analysis package for multi-speed gear box 
shafts- The program developed can also analyse shafts 
with circular cross section and general loadings. 

The salient feature of the program developed is, 
that a designer can quickly explore all the alternative 
design strategies by analysing them and select the best 
of all# which would otherwise be time csensuming. In 
case of gear box shafts, the designer can give different 
possible gear arrangements in gear box and t 

by analysing them through .iterations# can select a 
best feasible shaft. For this purpose the program 
offers several iteractive features lay which data can 
be modified and the entire design process can be 
repeated. 
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The purpose of shaft analysis is to determine the 
factor of safety of a shaft when given the geometry/ 
material/ and loading. In shaft design the two failure 
modes usually are yielding and fatigue/ which are governed 
by the deflectiop and stress along the shaft. The transfer 
matrix approach worhs effectively in finding the deflec- 
tion/ slope/ bending moment, and shear force at any 
location along the shaft. The differential equations 
used in the method are described in Chapter 2. Supports 
can be modelled as rigid or represented as sorings that 
may be translational or rotational in nature. Stress 
concentrations are added to more accurately model the 
shaft. However, stress concentrations are assumed to 
occur only at critical portion and not over the distri- 
buted area. Fatigue failure is considered at critical 
portions. The distortion energy or the maximum shear 
stress failure theory along with the Soderberg equation 
gives a reasonable prediction of failure. This analysis 
procedure is very general to enable modelling of shafts 
with complex loading and geometry. 

6.2 Recommendations for Further Worlc ; 

When required to design and analyse the gear 
box shafts, the designer has to give inputs to the 
program such as, the layout of the gears for each 
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arrangement. Keeping in mind the rules as to how to 
place them. These rules can be programed and thus one 
can get automatically the positions of the gears for 
each possible arrangement. Alternatively this layout 
of gear pairs and the location of support bearings can 
also be accomplished in an interactive graphical mode 
using a tablet or similar such input device. 

l-^Ihen the sizes of the shafts are finalized/ it 
is necessary to check that the operating speeds of the 
shafts do not coincide with the critical speeds of the 
shafts/ this part can be included in the presently 
developed program. 

The design of bearings at all support locations 
is also an important phase of the design process. Thus 
one can add in this vrork/ the necessary design calcula 
tions to select and size the appropriate type of bearing. 

The failure theories/ for brittle materials/ can 
also be added for specific analysis. 

Curve f i ttLng algorithms can be added for the 
stress concentration plots/ so that stress concentra- 
tion factors can be considered directly in a computer 
program. To get the optimized size of the shaft/ 
optimization methods can also be implemented. 
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APPENDIX-A 


DERIVATION OF TRANSFER MATRICES 


(A.I) Introduction; 


There are number of ways to obtain the transfer 
matrices from the equations of motion. The most direct 
of available methods is simply to rearrange the arbitrary 
constants of integration of solution in terms of the 
initial parameters [7,12]. Here only field matrix 
development is discussed. 


The fundamental equations of motion in first 


order form for the bending of a Timoshenko beam are. 


^ix 

Jtx 


© + 





(a. 1 ) 
(A. 2) 


= V + (K^-P) 9 - 


o 


% 

"■y 


c(Xr t ) 

(a. 3) 


= Ky + Q - W (x,t) (a. 4) 

s' 

These partial differential equations of motion for beams; 
in the case of steady-state motion, the loadings and 
responses (deflection, slope, bending moment, and shear 
force) all vary harmonically and reduces to ordinary 
differential equations as. 
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dy 

9x 

V 

= - e 

( A . D ) 


s 


dx 

M 

M 

“ El El 

(A. 6) 

dM 

dx 

= V + (K^ - P - r^.<*3^)© ~ 

c (A. 7) 

dx 

= (K - y - VJ 

(A. 8) 


These equations apply to the free vibration of a beam 
if the applied loadings C, W, are set equal to zero. 
Also, they are appropriate for static responses if c*. 7 is 
set equal to zero. 


(A, 2) Laplace Transformation: 

In matrix notation Equations (A. 5) to (A, 3) can 
be written as 


= AS + f 
dx 


and 


where. 



y 


0 


© 



S tr: 

M 

, f = 

-<=1^ lr> 


V 




0 -1 


IS. ^ 2 

0 K -P~^r tji> 


0 

l/EI 

0 



(A. 9) 


I/GA3 

0 

1 


. . (A- 10) 


0 


0 


0 
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The Laplace transform of Equation (A, 9) is written as 

s S(s) = S(0) + AS(s) -i- f (s) (a, 11) 

then 

S(s) = (IS - A)*"^ S(0) + (IS - A)~^ f (s) 

.... (A. 12) 

where I is a unit diagonal matrix. The inverse Laplace 
transform is 

S(x) = l“^ [(IS - A)“^]3(0) + "^[(IS - A)~%(s)] 

... (A. 13) 

The term L~^ [(IS - a)*"^] is the 4x4 transfer matrix, 
and the second term in equation (A. 13) is the fifth 
column of the extended transfer matrix. 

Now to find a field matrix for a Euler Bernoulli 
beam with shear deformation. Cnoose the case of a beam 
considering shear deformation with C, K, K , P, .jW/^1, 
r , equal to zero , taen 
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Combining equations (A. 13)/ (A.16)^. and (A«18) in the 
extended transfer matrix form 


-- 


1 


x^ 

X*" X 




y 


2 El 

6eT ^ GA 

s 

24EI 2"ga 

s 


y 

© 


0 

1 

X 

2 

X 



Q 


El 

2 El 

6 El 



= 





W.x^ 



M 


0 

0 

1 

X 

i 

~ 2 


M 

V 

1 

1 

j 

1 

j 

0 

0 

0 

1 

- W^x 


V 

1 

j 

X 

0 

0 

0 

0 

1 


1 


(A. 19) 


Similarly it is possible to derive the point matrix by 
this Laplace transformation. 



APPENDIX- B 


S>JAMPLE OF TRANSFER MATRIX METHOD 
(B,1) Introduction ; 

The transfer matrix method is capable of 
modelling complicated beams or shafts. An exarrple of 
this method's solution procedure# to a statically inde- 
terminate beam with an in-span rigid support is given 
in the following section below. The procedure for 
incorporating in-span rigid supports is complex, but 
straightforward in application, 

(B.2) Transfer Matrix Approach ; 

The purpose of this section is to give a detailed 
example with sequential steps for the transfer matrix 
method. The e>cample given in Fig. (B.l) is a cantilever 
beam with a rigid in-span support and a point load at 
the overhung end. The objective here is to determine 
the deflection, slope, bending moment and shear force 
at X = 0, X = Lj and x = L. The major objective is to 
determine the state vector at left end (x = O), because 
then the deflection, slope, bending moment, and shear 
force can be calculated at any point along the beam. 
First, the beam must be modelled in terms of sections 
that connect abrupt geometric changes or point 
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occurrences. Two sections exist for the e 2 :araple x = 0 
to X = L and from x = to x = L/ since a rigid in- span 
support occurs atx=L^ z as shown in Fig, (B.l). 
Assuming, no geometric changes occur along the beam and 
all the material properties are known, the transfer 
matrices along the beam can be formed. To apply the 
boundary conditions, the global matrix must be formed 
from the following equation, 

Oop [up = [u] 

(B.l) 

In this equation [u^] and [U 2 ] are field 
matrices and [u.] is point matrix for rigid in-span 
support at X = and [u^] is point matrix for the 
point load W. 

The field matrix [u ] is for x = o to x = 
and form eg. (2.15) is equal to 



(B.2) 
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Since shear deformation is vety small compared th the 
3 

term ~ L^/ 6 SI, thus is neglected in present example. 


The terms in equation (B,2) F , 

y ® 


, Fj^ and 


F^ are equal to zero sioce the beam is not loaded with 
uniformly distributed or linearly distributed loads. 
Thus the field matrix [u^] for x = o to x = L. is 


[Uj] = 


0 

0 

0 

0 


-L. 

1 

0 

0 

0 


~L^/2EI 

L,/EI 

1 

0 

0 


L^/2EI 


^1 

1 


0 

0 

0 

0 

1 


(B,3) 


Similarly the field matrix [u^] for x=L^tox=L can 
be given as 


[U3] = 


0 

0 

0 

0 


1 

0 

0 

0 


-L 2/2 El 

L2/EI 

1 

0 

0 


-L 2 / 6 EI 

L2/2EI 


1 

0 


0 

0 

0 

0 

1 


(B.4) 


[ U 2 ] is the point matrix for the concentrated load at 
the end of the beam and is given by using Table (2.2) 
as/ 
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1 0 

0 1 

0 0 

0 o 

0 0 


0 0 0 

0 0 0 

10 0 
0 1 -w 

0 0 1 


(B.5) 


To determine [ which is a point matrix for the 

rigid support, the subglobal matrix from x = to x = 
must be calculated. The subglobal matrix is egual to 
[ U 2 ] times [ U 2 ] and becomes 


[0j][u2] 


1 -L^ 

0 1 

0 0 
0 0 
0 0 


-L2/2EI 

L2/EI 

1 

0 

0 


--l|/6EI 

L2/2EI 

^2 

1 

0 


o 

o 

0 
-w 

1 


(B.6) 


L 


After obtaining the subglobal matrix, one of the boun- 
dary conditions at the right end is applied to calculate 

the point matrix [u. ] for in-span support in Table (2.3). 

^ * 

Since beam is free at the right end ' the boundary 
conditions at right end are that the bending moment and 
shear force are equal to zero (M^_l = 0 and = 0), 

Thus letting ^ i^ Ttable (2.3), [ U^] can be foxind 


to be as 
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Bv multiplying all the transfer matrices from left to 
right, the global matrix is equal to 

[u] = [ U2] 1U2] [ u^] [u^] 

which becomes as 



(B.9) 


Nox/, the remaining boundary condition at the right end 

(M =0) and the bounder'- condition at the rigid in- 
x=C» 

span support are applied as outlined in Section (2.S). 
The following equations are formed 
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\=L = 0=^Mn ^'o + ^14 


(B.ll) 


The terms, and P^ , come from the subglobal 

matrix from x = O to x =3 , and in this case this 

subglobal matrix is substituting the desired terms 

from the subglobal matrix, equation (B, 10) becomes 

0 = (- L^/2EI) + (- L^/ 6 EI) (B.12) 

The elements from the global matrix U,„,, U.,,, and P^, , 
are substituted into equation (B.ll), which leads to 


0 = Mo t Vo Li + L 2 W 


(B.13) 


After simultaneously solving equations (B,12) and (B,13) 

for the unknown state vectors, M and V become 

' o o 


- 3L^ M/2L^ 


(B. 14 ; 




(B.15) 


And since the boundary conditions at left end which is 
fixed, are Vo = 0 and = 0 j. the state vector at left 
end X o can be given by using above equations (B,14) 
and (B,15) and boundary conditions, as 


(4 


! W/2 

- 3L, W/2L^ 


(B.16) 



125 


The deflection/ slope/ bending moment/ and shear force 
are calculated at the locat'ons of interest by multi- 
plying the reouired matrices together v/ith the state 
vector at the left end. For example/ to find the def- 
lection, slope, bending moment, and shear force at just 
to the left of X = , the following equation is deve- 

loped. 



or 


y 


1 


-L“/2EI 

9 


0 

1 

L^/EI 

M 


0 

0 

1 

V ^ 

1 

j 

0 

0 

0 


1 

1 1 

ix=L. 

0 

0 

0 


S s 1 

j x=o 


(B,17) 

-L£/6EI 

0 


0 

L^/2EI 

0 


0 


0 


L 2 W /2 

1 

0 


-3L2W/2L^ 

0 

1 


1 


This gives 



-W L^L2/4EI 

- W 1^2 

-3L2 W/2L^ 


(B.IS) 


(B.19) 


0 
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The other location of interest is at x = L. The state 
vector at this point is found by using eauations/ 



[ 


[Uj] [ Uj] [Uj] 1 

(B.20) 

or 





% f x=L 

[ 

03 ^ 

[tTj] .'[ uj] [ 

(B.21) 


by using equation (B«21)^ the successive state vectors 
can be calculated just by using tlie previous state 
Vectors and proper transfer matrices in between them. 

Thus the state vector at x = L is given as 

W Li2/4EI + L 2 W/3EI 

-W h^/4El - L 2 W/2EI (B.22) 

0 
0 
1 

The in between state vectors can be computed by adjusting 
the length coordinates in the field matrix for the 
perticular sections and by properly arrangirg the 




transfer matrices 



